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Abstract 

We  introduce  a  new  mathematical  design  methodology  for  families 
of  codes,  suitable  for  different  applications  in  Radar  and  Communi¬ 
cations.  The  techniques  are  derived  from  Mulitiresolution  Harmonic 
Analysis  as  well  as  a  general  result  of  Benke  [2],  regarding  Rudin- 
Shapiro  Polynomials.  Some  explicit  estimates  for  the  performance  of 
the  algorithms  are  also  given. 


Motivation  and  Background 
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Since  the  classic  work  of  Golomb  [5],  numerous  mathematical  techniques 
have  been  devoleped  for  the  generation  of  code  families  with  high  auto-  and 
cross-correlation  performance.  With  few  exceptions  though  (e.g.  [3],  [1]), 
researchers  in  the  held  have  been  interested  in  the  developement  of  the  alge¬ 
braic  aspect  of  the  theory  of  code  design.  The  goal  of  the  present  paper  is  to 
demostrate  how  analytical  techniques  can  provide  higher  flexibility  to  build 
different  requirements  into  the  algorithms.  In  order  to  achive  this,  we  propose 
here  some  techniques  which  make  extensive  use  of  tools  from  Mulitiresolu- 
tion  Harmonic  Analysis  as  well  as  a  general  result  of  Benke  ([2])  regarding 
the  fundamental  construction  of  Rudin-Shapiro  Polynomials.  More  specifi¬ 
cally,  we  consider  codes  and  code-families  as  finite  approximants  of  bases  in 
inknite-dimentional  function  spaces  (e.g.  L2([0,1])  and  L2(M)).  An  approxi¬ 
mation  scheme  is  exhibited  with  the  desired  asymptotic  properties.  We  then 
obtain  a  variety  of  new  code  generation  algorithms  and  provide  some  explicit 
estimates  for  their  performance. 

The  paper  is  therefore  naturally  divided  in  three  parts:  a  first  part,  contain¬ 
ing  a  short  review  of  relevant  results,  the  second,  providing  statements  and 
proofs  of  the  main  theorems,  and  a  third  one,  with  numerical  results  a  short 
discussion  of  future  work. 


1  Algebraic  and  Anlaytic  Code  Design 

To  motivate  our  approach  to  the  design  problem  of  coding  sequences,  we  will 
follow  the  historical  development  of  the  subject  through  the  theory  of  Shift 
Register  Sequences.  Solomon  Golomb’s  classic  text  was  our  guide  in  this  [5]. 
Let’s  begin  with  some  preliminaries. 

1.1  Shift  Register  Sequences 

Linear  shift  registers  are  very  important  for  the  algebraic  theory  of  error- 
correcting  codes.  In  a  shift  register  (SR),  eventually,  a  sequence  will  repeat. 
This  is  because  for  a  binary  SR  sequence,  there  are  only  2r  possible  states 
(either  on  or  off,  for  each  tube).  So,  a  repetition  occurs  in  the  first  2r  states. 
However,  we  can  improve  on  that  bound,  since  if  we  have  a  state  of  all  0’s, 
the  shift  register  will  continue  producing  0’s,  which  means  its  period  is  just 
1.  So,  the  period  of  a  binary  shift  register  is  at  most  2r  —  1. 
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Definition  A  sequence  generated  by  an  r-tube  shift  register  will  be  said  to 
have  maximum  length  if  its  period  is  p  =  2r  —  1. 

Lemma  1.1  Any  r  inputs  and  r  outputs  of  a  maximum  length  r-tube  shift 
register  sequence  completely  determine  all  of  the  outputs. 

Proof  Consider  r  inputs  {ai,  a2, . . . ,  ar }  and  r  outputs  {bi,  b2, . . . ,  br}. 
Then, 


Ip  =  Ciar  ©  c2ar_i  ©  •  •  •  ©  crai  (mod2) 
b-2  =  c\b\  ©  c2ar  ©  •  •  •  ©  cra2  ( mod2 ) 


br  =  Ci6r_i  ©  c26r_2  ©  •  •  •  ©  crar  (mod2) 


So,  we  want  to  find  the  recursion  coefficients,  {ci,  c2, . . . ,  cr}.  Hence,  in 
matrix  form  this  system  becomes: 


Hence, 


/  h  \ 

b2 

= 

\  br  ) 

CLrp  CLy _ | 

r 


b\  a, 

V  br- 1  h-2 


&\ 


\ 


) 


(  Cl  \ 

C2 


\Cr  ) 


/c\ 

(  ar  ar- 1  ■  ■  •  a\  ^ 

- i 

f  h\ 

c2 

= 

b\  ar  ■  ■  •  o2 

b'2 

\  Cr  ) 

\  br-i  br-2  ■■■  ar  J 

\br  ) 

if  the  inverse  exists. 


Assume  that  it  does  not.  Then,  there  is  a  row  such  that  it  is  a  Z2  multiple 
of  another  row.  Hence,  it  is  either  equal  to  that  row  or  is  a  row  of  0’s. 

Say  it  is  the  latter.  So,  say  the  fc-th  row  is  all  0’s.  If  k  ^  1,  then 

b\  *  bfc—  i  bfc  ar  ■  Ufc  0. 


But,  this  means  that 
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bk+i  —  c\bk  ©  •  •  •  ©  Ckb\  ©  Ck+iar  ©  •  •  •  ©  crak+i  —  Oci  ©  •  •  •  ©  0 cr  —  0. 

In  the  same  manner,  we  see  that  all  the  bk  are  0,  meaning  that  we  have 
a  blank  tape,  contradicting  that  the  shift  register  has  maximum  length.  If 
k  —  1,  then  all  the  a*  are  0,  meaning  that  all  the  inputs  are  0.  This  implies 
that  all  the  outputs  are  0,  which  is  again  a  contradiction. 

Now,  say  that  there  is  a  row  equal  to  another  row.  Let  the  j-th  row  equal 
the  k  +  j-tfi  row.  Then,  it  can  be  seen  that  the  elements  of  the  sequence: 
{cq,  a2, . . .  ,ar,bi,b2,  ■  ■  ■ ,  br}  repeat  every  k  elements.  Hence,  the  sequence 
has  period  k,  contradicting  that  it  is  of  maximum  length,  since  k  <  r  <  2r— 1. 


So,  the  matrix  is  invertible.  Hence,  we  can  solve  for  ci, . . . ,  cr.  Now,  construct 
the  matrix: 


/  ci  1  0  •  •  •  0  \ 

U  o  1  0 


\Cr  0  0  •  •  •  1  ) 

Then, 


(  Cl 

1 

0  •• 

•  0  \ 

\ 

C2 

0 

1  •• 

•  0 

Q"nQJn-\- 1  *  *  *  ^n+r  ) 

\  Cr 

0 

0  •• 

•  !/ 

(  ®n+r+l®n+r+2  '  '  '  ®n+2 r  ) 


for  all  n.  Hence,  this  matrix  yields  the  entire  output.  § 


Remark  To  see  that  the  general  case  (for  Zp)  is  also  true,  we  use  the  same 
argument,  only  this  time,  we  assume  that  some  row  is  a  Zp  multiple  of 
another.  So,  say  that  Vj  is  the  j-tfi  and  vk+j  is  the  k  +  j-th  row.  Then, 
Vj  =  dvk+j  (mod  p).  So,  we  see  that  after  every  k  outputs,  the  next  k  outputs 
are  d  times  the  previous  set  of  k  outputs.  Now,  by  Fermat’s  Little  Theorem, 
dP~l  =  1  (mod  p).  This  means  that  after  k(p  —  1)  places,  the  sequence 
repeats.  But,  k(p  —  1)  <  pr  —  1.  Hence,  the  sequence  is  not  of  maximum 
length.  This  is  a  contradiction.  Hence,  the  rows  are  linearly  independent 
and  the  rest  of  the  argument  from  above  follows. 
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In  a  variety  of  engineering  applications,  there  arises  a  need  to  have  ”  random” 
sequences.  Since  a  computer  is  a  finite  state  machine,  true  randomness  can 
not  be  produced  on  it.  Hence,  there  is  a  need  to  produce  sequences  that 
appear  random.  A  good  model  for  binary  random  sequences,  is  flipping  a 
fair  coin.  From  statistics,  there  are  certain  things  one  would  expect  from 
such  a  model: 


•  The  number  of  +l’s  (heads)  is  about  the  same  as  the  number  of  —  l’s 
(tails). 

•  Short  runs  (consecutive  streaks  of  heads  or  tails)  are  more  likely  to 
occur  than  long  runs.  Precisely,  half  the  runs  have  length  1,  one  fourth 
have  length  2,  one  eighth  have  length  3,  etc. 

•  There  is  also  a  certain  property  about  the  autocorrelation  of  such  se¬ 
quences.  Autocorrelation  measures  how  similar  a  sequence  is  to  a  shift 
of  itself.  One  would  expect  that  the  autocorrelation  peaks  at  no  shift 
(being  identical  to  itself),  and  is  smaller  for  positive  shifts. 

The  algebraic  definition  of  a  random  sequence,  which  in  engineering  is  called 
noise,  is  then: 

Definition  We  say  that  a  binary  shift  register  sequence,  A  =  {oi,  a,2, . . .}  of 
period  p  is  a  pseudo-noise  sequence  if: 

1. 

p 

\yan\<i 

n=  1 


2. 

3. 


For  every  run  of  length  k ,  there  are  two  runs  of  length  k  —  1 
The  cyclic  autocorrelation  function  a  has  the  property: 

p  if  r  =  0 

n=l 


pa(A)(r)  =  Y  anan+T 


O  if  0  <  r  <  p 


The  three  randomness  postulates  all  reflect  the  model  of  randomness  de¬ 
scribed  above.  The  third  postulate  says  that  the  cyclic  autocorrelation  func¬ 
tion  is  two  valued.  This  comes  from  the  fact  that  all  shift  register  sequences 
are  periodic.  It  turns  out  that  all  maximum-length  shift  registers  satisfy  the 
three  postulates.  The  proof  of  this  fact  is  not  hard,  but  is  beyond  the  scope 
of  the  present  discussion  and  may  be  found  in  Golomb’s  text. 
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2  Crosscorrelation  and  Autocorrelation  Prop¬ 
erties 

After  the  introductory  study  of  SR  sequences,  we  proceeded  to  look  at  Dilip 
Sarwate’s  1979  paper  Crosscorrelation  and  Autocorrelation  of  Sequences  [6]. 
This  provides  a  good  feel  for  some  basic  analytical  aspects  of  the  project. 

Let  A"  be  a  family  of  K  complex  sequences  of  period  N.  For  sequences, 
u,v  G  X,  the  periodic  crosscorrelation  function  7 (u,v)(-)  is  defined  as: 

n-  1 

7 (u,v)(l)  =  Y,  uiyi+i 

i= 0 

and  the  periodic  autocorrelation  function  «(«)(•)  is  dehned  as: 

N- 1 

a(M)(0  =  y  UiUi+i 

i= 0 

We  now  prove  a  technical  lemma: 

Lemma  2.1  Let  u,v  be  two  complex  sequences  of  length  N.  Then, 

N-l  N- 1 

y  |7 (u,v)(l)\2  =  ya(u)(l)a(v)(l) 

1= 0  z=o 
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Proof 


TV-1 


TV— 1 


E^ItOwXOT  =  2^y(u,v)(Z)'y(u,v)(Z) 


1=0 


Y1  uivi+l 


1=0 

TV— 1  /  /TV-1  \  / N-l 

=  E  E Ui 

«=o  y\i=o  /  \j=o 

TV-1  TV— 1  TV-1 

=  ESE  V’iVi-^-iUjVj+i 

1=0  i= 0  j=0 
N-l  TV-1  TV-1 

=  EEE  VtjVi^i  tlj  VJ+1 

i= 0  j=0  T=0 
TV-1  TV-1  TV-1 

=  EE^E 

i=0  j=0  Z=0 

Since  the  sequences  are  periodic  with  period  N,  we  can  work  modulo  N.  So, 
letting  rn  —  j  —  i  and  k  =  i  +  /,  we  notice  that: 


TV-1 


ItOwxop 


1=0 


N-l  N-l  N-l 

Vi+lVj+l 

i=0  j= 0  T=0 

TV— 1  TV— 1  TV-1 

EE  E  t’Wfc+m 

i=0  m=0  fc=0 

W  ( (s2  — 

/  ^  I  I  /  j  Uji  U>i+m,  I  I  /  ^  'VkVk+m 

m= 0  y  \  i=0  /  \  k= 0 

TV— 1 

a(u)(m)a(v)(m).  § 

m= 0 


We  are  now  ready  to  state  the  main  result  of  Sarwate’s  paper.  Dehne  the 
following  two  numbers  A  and  T : 

T  =  max{|7(rz, v)(7)|  :  u,v  G  X,  0  <  l  <  N  —  1} 

A  =  max  {\a(u)(l)\  :  u  G  X,  0  <  l  <  N  —  1} 


Theorem  2.2  (Sarwate,  1979)  For  any  set  X  of  K  sequences  satisfying 
a(w)(0)  =  N  for  all  u  G  X,  the  following  holds: 


/r2\  N  -  1  fA2\  ^ 
ViV  y  +  N(K  -  1)  \  N  J-1' 


Proof  For  all  u,  v  G  X, 


N-l  N—l 

E  E  E  i7(«,f)(or  +  E  E  i“(«)(oi2 

u£X  v£X  1=0  u£X  1=0 

u^v 


N-l 


1  t(«,u)(z)|2 

u&X  v£X  1=0 
N-l 

EEE  a(u)(l)a(v)(l)  (by  Lemma  2.1) 


ueX  v£X  1=0 
N-l 


J^a(u)(0  j  J^a(v)(o 

i=o  y  \uex  J  \vex 

N-l 

EiE“(“)<i)i2 

=0  u€X 

\  2  N—l 

E  “(“)(<>)  +EiE“(“)(')ii 

uex  /  i=i  uex 

N—l 

a-2jv2  +  E  lE“(“)WI2 


1=1  u&X 


since  there  are  K  elements  in  X  and  for  each,  a(tt)(0)  =  N  by  assumption. 
Now,  for  the  first  term  on  the  left  hand  side,  there  are  K(K—  1)  different  pairs 
of  sequences  in  X.  For  the  second  term,  for  l  =  0,  there  are  K  sequences, 
with  a(tt)(0)  =  N  for  each  u  G  X.  Hence,  the  LHS  is  bounded  from  above  by 
K ( K  —  1) AT2  +  K N2  +  K (N  —  1)H2.  On  the  other  hand,  the  right  hand  side 
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is  bounded  from  below  by  K2N 2  since  the  second  term  is  always  positive. 
Hence, 


K(K  -  1)AT2  +  KN2  +  K(N  -  1  )A2 
=>  K(K  -  1)AT2  +  K(N  -  1  )A2 
/r2\  (n  —  i)  (A2\ 
\NJ  N{K  -  1)  V  N  ) 


>  k2n 2 

>  KN2(K  —  1) 

>  1-  I 


Definition  We  say  that  a  code  is  a  complex  sequence  where  each  element 
has  norm  1. 

Example  Consider  the  family  of  codes  of  prime  length  L:  U  =  {u^\u^2\  . . 
where  for  all  N ,  }fc=o-  So,  let’s  calculate  the  constants  T  and 

A. 


a(u(N))(l )  = 

L- 1 

V — >  27 rikN  27t  i{k+l)N 

L  e  L 

k= 0 

L- 1 

\j — >  2nikN  2ni(k+l)N 

)  ,(  '■  < 

— 

k= 0 

L- 1 

y — -\  2-kHN 

= 

k= 0 

2t rilN 

Le  L 

Hence,  A  =  L.  As  for  T, 

7 

7,-1 

v — >  27 rikN  27ri(fc+Z)M 

=  l  e  l 

k= 0 

L— 1 

V — >  2-KikN  2irt(fc  +  i)M 

=  l^e  L  e  L 

k= 0 

L—l 

^  2iri(kN  —  kM  —  IM) 

=  Le 

k= 0 
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Since  the  codes  are  periodic  with  prime  period  L,  we  can  work  Modulo  L. 
Hence,  as  we  vary  k,  kN  —  kM  —  IM  eventually  cycles  through  all  of  Z l. 
Then,  this  sum  is  just  the  sum  of  the  Lth  roots  of  unity.  Hence,  they  sat¬ 
isfy  the  equation  ujl  =  1.  From  basic  algebra,  this  sum  is  —  where 

aL-i  and  cll  are  the  coefficients  of  the  (L  —  l)th  and  Lth  terms  of  the  poly- 

_ — ,  r  _1  2-KiikN  —  kM  —  IM) 

nornial.  Therefore,  Yhk= o  e  L  =  0.  Hence,  T  =  0. 

Lemma  2.3  Let  X  be  a  family  of  K  codes  of  period  N,  as  in  Theorem  2.4, 
s.t.  A  =  N  and  T  =  0.  Then,  K  <  N 

Proof  Let  A  —  N,  T  =  0.  Then,  by  Sarwate’s  Theorem,  >  1-  Hence, 
N  >  K.  | 

Lemma  2.4  If  a  family  of  K  codes,  X,  has  the  property  that  T  =  0,  then, 
there  exists  a  code  u,  such  that  a(u)(l)  =  N  for  some  l,  0  <  l  <  L  —  1. 

Proof  By  Sarwate’s  Theorem,  if  T  =  0,  A2  >  N  Now,  if  we  represent 

our  codes  as  vectors  Vj,  then  we  note  that  if  T  =  0,  for  any  Uj ,  14.  (j  7^  k ), 
(uj  ,  U]f)  =  0.  Hence,  all  of  the  vectors  are  linearly  independent.  Therefore, 
K  <  N .  Thus,  we  see  that  A2  >  N 2 .  Hence,  A  >  N.  But,  we  also  know  that 
A  <  N  since  the  maximum  autocorrelation  occurs  at  zero-shift.  Therefore, 
A  =  N.  Thus,  3 u  G  X,  l  >  0,  s.t.  a(u)(l)  =  N  .  | 

3  Fourier  Transforms 

Fourier  Transforms  serve  an  important  role  in  signal  processing.  A  code 
can  be  considered  as  a  discretely  sampled  version  of  a  signal  of  constant 
amplitude.  So,  we  can  extend  the  use  of  Fourier  Transforms  to  codes.  So, 
for  discrete  codes,  we  will  make  use  of  Discrete  Fourier  Transforms: 

Definition  Consider  a  complex  periodic  code  of  period  N\  {rco,  X\, . . . , 

We  define  the  Discrete  Fourier  Transform  as  the  complex  periodic  code  of 
period  N:  {y0,  X 1 ,  -  -  -  ,  Xjv-i},  where  for  each  n  —  0, 1, , . .  N  —  1, 

1  N -1 

A  ^  "\  —27 rikn 

=  Xke 
k= 0 
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Proposition  3.1  Let  C  =  {c0,  C\% . . . ,  cjv-i}  be  a  complex  periodic  code  of 
prime  length  N .  Let  C  =  {co,  ci, . . . ,  cat_i}  be  the  DFT  of  the  code.  Then, 

a(C)(l )  =  0 


for  all  l  =  1, . . . ,  N  —  1. 

Proof  For  each  n  =  0, 1, . . . ,  N  —  1,  cn  =  ^2k=o  Cke  n  .  Hence, 

N—l 

C^(C)(^)  ^  ^  CmCm+i 

m= 0 

JV-l  /  /JV-1  \  /jv^l  \ 

1  X — ■>  /  /  ^ —2-TTikm  \  I  \ — A  —2irij(m+l)  | 

=  Z>e  "  2>e^^ 

m= 0  y  \k=0  /  yi=0  J 

,  jV-lV-liV-l  _ 

—  2-nikm  —  2n 

Cke  N  Cje  n 

k= 0  j= 0  m=0 

iV-lV-l  V-l 

1  V  ^  \  ■*,  _  \  -V  2n i(jm+jl  —  km) 

=  Jp2^z^  CkCi  e  N 

k= 0  j= 0  m= 0 

Note  that  for  each  j,  k  the  third  summation  is  just  the  sum  of  the  IVth  roots  of 
unity  as  before  since  N  is  prime.  Hence,  it  is  0.  Therefore,  a(C)(l)  =  0.  j§ 


We  can  extend  the  theory  of  discrete  codes  to  the  theory  of  continuous  ones 
(ie:  L2(M),  Cp ,  etc.).  Here,  we  will  make  use  of  the  continuous  version  of 
the  Fourier  Transform  extensively.  Hence,  we  proceed  to  extend  the  discrete 
results  from  Section  2  to  continuous  analogues. 

The  identity  in  Lemma  2.1  is  true  in  the  discrete  case.  But  what  happens 
in  the  continuous  case?  We  have  a  natural  way  of  defining  crosscorrelations 
and  autocorrelations  for  continuous  signals  (or  functions),  namely  via  the 
convolution  product.  So,  if  /,  g  G  L2(M),  we  define  the  crosscorrelation  of 
two  functions  as: 

r+oo  _ 

(f*9)(y)  =  (f*9)(y)=  f(x)g(y  +  x)dx 


and  the  autocorrelation  as: 


_  r+oo  _ 

(/  *9)  {y)  =  {f*f){y)  =  /  f(x)f(y  +  x)dx, 

where  for  any  function  /,  f(x)  is  defined  as  f(x)  =  f(—x). 

We  also  define  the  Fourier  Transform  of  a  function  as: 

_  p+00 

/K)  =  ^[/W]=  /  f(x)e-2’^dx. 

J  —OO 

So,  is  it  true  that  (/*/  ,  g*  g)  =  ||/*g||2?  The  answer  is  yes. 
Lemma  3.2  Assume  that  f,g  are  complex  functions  in  L2(M).  Then, 

(/*/  .,  g*g)  =  ||/*3||2, 

where  the  inner  products  and  norms  are  taken  w.r.t.  L2(M). 

Proof  The  Fourier  Transform  is  an  isometry.  So, 

\\f*g\\ 2  =  (f*9,f*9) 

—  (fg  ,  fg)  (by  the  Convolution  Theorem) 


First,  we  note  that, 


/«)  =  if(x)  ,  e2*“£} 
=  (e-2*“£  ,  /(*)} 


=  (/(x)  ,  e  2ma;?) 


r»+oo 


f(x)e2nix ^  dn(x) 


r»+oo 


f(—y)e2nd-y)^  dy(y)  (by  changing  variables  x  =  — 
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Hence,  /(£)  =  /(£)■ 


ll/^ll2  = 


r+o o 

"  — oo 
/*+oo 


f  (Of  (ZMZMt)  MO 


=  {if ,  99) 

=  (fi ,  <?f) 

=  (f  *1 ,  9*9) 

=  (f  *  f  1  9*9)  (since  the  F.  T.  is  an  isometry). 


I 


Furthermore,  assume  that  they 

Now  suppose  instead  of  a  function  in  L2(M),  we  have  a  function  in  L2([0, 1)). 
Can  we  prove  a  comparable  result?  The  answer  is  yes.  All  of  the  same 
arguments  apply  as  in  the  case  of  L2(M).  The  only  thing  we  have  to  modify 
is  that  we  note  that  the  map  /  1— >•  /  (where  /  is  the  Fourier  expansion  of  /)  is 
an  isometric  isomorphism  between  L2([ 0, 1))  and  Z2( Z)  by  the  Riesz-Fischer 
Theorem.  Then,  we  only  require  that  the  Fourier  series  of  /  and  g  converge 
uniformly  so  that  we  can  interchange  the  infinite  sum  and  the  integral,  and 
all  of  the  previous  arguments  apply. 


Continuing  in  extending  our  results  from  the  previous  section  to  spaces  of 
non-discrete  functions,  such  as  L2,  we  now  ask  if  we  can  establish  a  bound 
similar  to  the  Sarwate  one  above.  Define  the  following  two  numbers  A  and 

f: 


A  =  max  {Hu  *  u||sup, 

T  =  max  {||w  *  u||SUp,  u,v  E  X,  u  ^  v} 

The  next  Theorem  gives  such  a  result. 


Theorem  3.3  Consider  a  family  X  C  L2([0, 1))  of  K  continuous  functions 
s.t.  for  all  u  G  X,  \  J( u(x )  dfi(x)  \  =  1.  Then, 


A2 

K 


+  (K  -  1) 


f2 

K 


>  1. 
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Proof 


EE  ||u  *  v\\2  +  ^  \\u  *  u\ 


ueX  veX 
u^v 


ueX 


\Ukv\ 


EE 

uex  vex 

EE  (ukU  ,  Vkv) 

uex  vex 

UkU  ,  V  k  v) 

uex 

II E 


vex 


U  k  U\ 


uex 


Now,  as  in  Lemma  3.2,  we  note  that, 


(u 


k  u)(x)  =  X\\u(k)\2e2'Rlkx ,  and  u(0)  —  /  u(x)  d/i(x ) 

J  o 


fcez 


Hence, 


E 


2iir'2 


«ew  vex 


UkU  ||  = 

uex  uex  vex  kez 

=  K^)i2i^)i2  + 

u^X  v€.X  Ai£Z\0 

=  +  Kfc)i2i*Kfc)i2- 

uex  vex  kez\ o 


The  left-hand  side  is  bounded  from  above  by  K A2  +  K(K  —  1)T2  since 


l/(*)l  <  ll/ll 


sup 


\f(x)\2  dfi(x)  <  \\f(x)\ 


sup* 


On  the  other  hand,  the  right  hand  side  is  bounded  from  below  by  K 2  since 
the  second  term  is  always  non-negative.  Hence, 


K A2  +  K(K  -  l)f 2  >  K2 


A2  f2 

-  +  {K-  1)— >1.  I 
K  y  J  K  ~  " 
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Remark  We  can  easily  make  this  into  a  more  general  result,  where  instead 
of  continuous  functions  we  use  any  L2([ 0, 1))  functions  and  substitute  the  L2 
norm  for  the  sup  norm. 


4  Rudin-Shapiro  Sequences 

An  area  of  coding  theory  where  analytic  methods  are  used  extensively  is 
the  study  of  codes  generated  by  Generalized  Rudin-Shapiro  Systems.  In  our 
study  of  these  systems,  we  used  George  Benke’s  1994  paper  of  the  same  name 
[2]-  ‘ 

The  classical  Rudin-Shapiro  polynomials  are  trigonometric  polynomials  de¬ 
fined  recursively  as: 


Po(x)  =  1,  Qo{x)  =  1 
Pn+i(x)  =  Pn(x)  +  e2”i2nxQn(x ) 

Qn+i(x)  =  Pn(x)-e2”i2nxQn(x). 

Lemma  4.1  Let  Pn(x)  and  Qn(x)  be  as  above.  Then,  for  any  n  >  0,  we 
have: 

\Pn+i{x)\2  +  \Qn+i(x)\2  =  2(|Pn(a;)|2  +  \Qn(x)\2). 


Proof 


Pn+2{x)  | 2  +  \Qn+2(x)\2 


+ 


Pn+i(x)  +  e2”i2n+lxQn+l(x)  +  Pn+1(x)  -  e2-i2ri+la 
2(\Pn+l(x)\  +  \Qn+i(x)\  ) 


Qn+1  {x) 


e27r  i2n+lxPn+1(x)Qn+1(x)  +  e~2™2n+lxQn+l(x)Pn+i(x ) 
e2”i2n+1xQn+1(x)P^(x)  -  e-2"2"+1^Pn+iW 
^{\Pn+l{x)Y  -\-\Qn+l{x)\2).  | 


Lemma  4.2  For  any  n  >  0,  we  have: 

\Pn(x)\  <  V2VN,  and  \Qn(x)\  <  V2\/N 
where  N  =  2n  =  \\Pn(x)\\l  =  \\Qn(x)\\l. 
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Proof 

\Pn(x)\2 +  \Qn(x)\2  =  2(\Pn_l(x)\2  +  |Qn_i(a;)|2)  (by  4.1) 

=  2  *  2n_1(|P0(^)|2  +  |Qo(^)|2)  (applying  4.1  repeatedly) 
But,  |P0(a;)|2  +  |Qo(^)|2  =  2  Hence, 

\Pn{x)\2  +  \Qn(x)\2  =  2  *  2n. 

This  means  that  \Pn{x)  <  s/2y/N .  The  same  holds  true  for  Qn(x).  | 


XJx)  = 


U  = 


DJx )  = 


Alternatively,  we  can  write  this  system  in  matrix-vector  form.  Let 

Pnix) 

Qn  {x ) 

Hence,  Equations  (1),  (2),  and  (3)  become 

Xn+1  =  V2UDn(x)Xn(x). 


(  l/y/2  l/y/2  \ 

v  i  /  V-  —  i  /  V- ) 


0 

2iri2nx 


This  motivates  the  generalization  of  the  construction.  We  begin  to  generalize 
it  by  letting  e  =  {en},  n  =  1,2, ... ,  be  an  arbitrary  sequence  of  0’s  and  l’s, 

and  letting  F  =  ^  ^  q  ^  •  Then,  we  define  the  sequence  Xe>n+1  as: 
Xe.n+1(x)  =  V2UDn{x)FenX^n{x). 


So,  if  we  consider  all  the  possible  choices  for  e,  at  each  stage  n,  we  obtain 
2n  polynomials  with  2"  coefficients.  So,  we  form  a  matrix,  where  for  each 
polynomial,  there  is  a  row  of  its  coefficients.  Hence,  the  first  three  matrices 
of  the  construction  are: 


1  1-1 
1  1-1 
1-1  1 
1-1  1 
-111 
-111 
-1  -1  -1 

-1  -1  -1 


(l  1  1  -1\ 

11-11 
1-1  1  1  ’ 
V  i  -i  -i  -i  / 

1  1-1  1  \ 

-1  -1  1  -1 
111-1 
-1  -1  -1  1 
1  -1  -1  -1 
-1111 
1-111 
-1  1  -1  -1 


and 


,  respectively. 
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It  can  be  clearly  seen  that  the  rows  are  orthogonal.  The  first  two  rows  are 
the  Nth  classical  Rudin-Shapiro  coefficients.  We  shall  call  this  construction 
the  ”  append”  rule.  It  is  so  called  since  at  each  stage  of  the  iteration,  new 
higher-order  terms  are  added  to  the  end  of  the  polynomial.  A  related,  but 
different,  construction  is  given  by: 

Xn+i(x)  =  V2UT(x)Xn(2x). 

where 

T(X)  =  f  g  ^ 

This  is  called  the  interlace  rule,  since  new  terms  are  not  added  to  the  end, 
but  rather  ”  interlaced”  in  between  the  previous  ones.  As  with  the  ”  append” 
construction,  we  can  also  consider  the  entire  set  of  coefficient  blocks  generated 
by  all  possible  e.  As  in,  for  each  e, 

Xe,n+i(x)  =  V2UT(x)Fe”Xe,n(2x). 

Now,  we  can  generalize  these  constructions  to  arbitrary  pxp  matrices.  Notice 
that  for  the  ’’append”  construction,  for  n  =  1,  the  generating  matrix  was 

^ .  Now,  we  can  have  the  system  of  polynomials  be  generated  by 

any  pxp  matrix,  where  p  is  prime. 


So,  how  do  we  get  at  this  generalization?  Note  that  in  the  classical  Rudin- 

Shapiro  Append  rule,  given  the  original  matrix  (  °00  Qo1  ) ,  and  two  con- 

\  ®io  °n  ) 

secutive  rows  R% ,  Rj ,  we  form  four  new  rows  for  the  next  iteration: 


(ooo R-i  «oi Rj)  i  (oio R-i  On Rj)  i 


( (/no  / a§\Ri) ,  and  (cpo-ht;  On Ri)  ■ 

Note  how  the  rows  of  the  previous  matrix  are  cyclically  permuted  in  the  con¬ 
struction,  and  the  entries  of  the  original  matrix  appear  as  coefficients  in  the 
new  one.  Thus,  generalizing  this  construction  to  the  pxp  case,  we  see  that 
expressing  the  row  and  column  indeces  in  a  base-p  representation  vector  is 
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most  convenient.  Hence,  letting  u,  u  G  Z”  and  p  —  {l,p, . . . , pn~ 1 } ,  the  gen¬ 
eralized  Append  and  Interlace  polynomials  (AnjV(x)  and  IH:U(x),  respectively) 
are  given  by: 


An,i/  (%) 
In,v  (*^) 


where  all  subscript  addition  is  computed  modulo  p. 


These  generalized  Rudin-Shapiro  systems  have  interesting  applications  to 
coding  theory.  One  could  construct  a  code  using  the  coefficients  of  the 
trigonometric  polynomial.  How  good  are  these  codes,  especially  in  terms 
of  the  Sarwate  bound  proven  earlier?  First,  we  change  notation  slightly  and 
say  that  z7  G  is  now  a  code  given  by  the  coefficients  of  the  polynomial 

Lemma  4.3  Let  n  >  1  and  p  >  2.  Let  A{n )  =  ^  ,  T(n)  —  ^  at  the  nth 
stage  of  the  process.  Also,  say  that  A  =  \ffdU  for  some  unitary  matrix  U. 
Then,  for  each  n, 

A(n)  <  (/ 3/p)n  and  r(n)  <  (/3/p)n. 

Proof  For  each  l  >  1  ,a(u)(l)  <  a(tt)(0)  since  the  autocorrelation  function 
peaks  at  zero-shift.  Now  note  that  for  each  code  z7,  o(z7)(0)  =  \An^(x)\2  dx  = 

/ 3n  by  Theorem  5  in  [2],  Hence, 

A  max{a(P)(0),PgZg"}  / 

Pn  ~  pn  \p/ 

Now,  for  T,  observe  that  if  we  represent  our  sequences,  z7,  V  as  vectors,  we 
can  write 

7(z7,z/)(/)  =  (V \SiV), 

where  Si  is  the  linear  transformation  that  shifts  a  vector  cyclically  by  l  places. 
Since  Si  G  Sym(pn),  it  is  an  isometry  of  Mp”.  Hence, 

<  \W\\  ll^^ll  =  N|  ll^ll  =  PnA  IIAmHI2 
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But,  by  Theorem  5  in  [2], 

||Ai,i7|l2  IIAi,Hl2  =  Pn- 

Hence,  T(n)  =  ^  <  (J)  .  | 


Lemma  4.4  Let  A  =  /~j3U,  for  some  unitary  matrix  U  and  constant  (3.  Let 
V  be  the  code  given  by  the  coefficients  of  the  Append  polynomial,  An^(x). 
Then,  the  Discrete  Fourier  Transform  of  u  is  given  by: 


5  =  {^-(0),  LjkaA),  ....  ^An,A—f)} 


Proof  By  the  definition  of  the  DFT,  if  z7  =  {uq,  u\,  . , then 

v  —  {Xo,  Xi,  •  •  • ,  Xpn~i}i 


where  for  each  m, 


Xm 


pn- 1 

—  l-nihra 

uke  pn 

k= 0 


But,  note  that  for  each  m  this  sum  is  just  the  Append  polynomial  evaluated 
at  x  =  —m/pn  since  the  coefficients  uk  are  just  the  coefficients  of  Anp(x) 
taken  in  the  same  order.  Hence, 


"  pn^nA  pn  )>  •  •  •  >  pn^n,p(  ^ 


P 


)}•  I 


Lemma  4.5  Let  u  be  as  above.  Then,  for  each  m  —  0, 1, . . .  ,pn  —  1, 

Proof  By  Corollary  4  in  [2],  for  all  x ,  \Anj(x)\  <  s/p  (vy?)”  •  Substituting 
into  the  above,  we  see  that  for  all  m, 

lx™l  =  <  p^p{F)‘  =  Vp(jP) 


I 
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Theorem  4.6  Say  that  f  is  a  function  in  L2((0, 1]),  such  that  for  all  k, 
3 C  >  0,a  >  1  such  that 

i/wi  <  F, 

an 

where  f(k)  is  the  k-th  Fourier  coefficient.  Furthermore,  assume  that  the 
Fourier  series  has  no  more  than  pn  terms,  where  p  is  some  prime.  Then, 
||/*/||i  ->  0  as  n  ->  oo. 

Proof  Say  the  Fourier  series  of  /  has  pn  terms.  /  *  /  =  Ylk=i  \f(k)\2e2mkx , 
as  before.  Then,  we  have  that 

pn 

f*f  <  sup  |/(fc)|2Ve2^ 

k<P-  k=1 

C 

<  Dvn(x ), 

-  an  P  V  J’ 

where  Dpn(x )  is  the  Dirichlet  kernel.  From  standard  Fourier  Analysis,  we 
know  that  ||-Dp™(:r)||i  =  4^  logpn  +  0(1)  (consult  [4]  for  a  complete  proof). 
Hence, 

II/* /Hi  <  ^(nlogp  +  0(l)) 

an 

-»•  0  +  0  =  0 


as  n 


(X). 
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The  above  theorem  has  important  consequences  for  us.  Note  that  for  any 
code  v  that  satisfies  the  hypotheses  of  4.4,  we  can  see  that  for  each  k,  the 
k-th  Fourier  coefficient  of  the  periodic  autocorrelation  function  a(z?)  is  given 
by  |  Xfc  | 2  where  Xk  is  the  k-th  Fourier  coefficient  of  our  code.  Now,  taking 
the  Inverse  DFT  of  the  DFT  of  the  autocorrelation  function,  we  can  recover 
the  autocorrelation  function  itself.  But,  notice  that  we  can  use  an  argument 
similar  to  the  one  in  the  preceding  theorem  to  get  a  bound  on  the  IDFT. 
Hence,  we  see  that  with  binary  codes  given  by  the  classical  RS  construction 
as  n  — >  cx),  a(y)(l)  — >  0  for  all  /  <  2n. 
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